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Linear optics underpins tests of fundamental quantum mechanics and computer science, as well 
as quantum technologies. Here we experimentally demonstrate the longstanding goal of a single 
reprogrammable optical circuit that is sufficient to implement all possible linear optical protocols 
up to the size of that circuit. Our six-mode universal system consists of a cascade of 15 Mach- 
Zehnder interferometers with 30 thermo-optic phase shifters integrated into a single photonic chip 
that is electrically and optically interfaced for arbitrary setting of all phase shifters, input of up to 
six photons and their measurement with a 12 single-photon detector system. We programmed this 
system to implement heralded quantum logic and entangling gates, boson sampling with verification 
tests, and six-dimensional complex Hadamards. We implemented 100 Haar random unitaries with 
average fidelity 0.999 ± 0.001. Our system is capable of switching between these and any other 
linear optical protocol in seconds. These results point the way to applications across fundamental 
science and quantum technologies. 


Photonics has been crucial in establishing the foun¬ 
dations of quantum mechanics [T], and more recently 
has pushed the vanguard of efforts in understanding new 
non-classical computational possibilities. Typical pro¬ 
tocols involve nonlinear operations, such as the genera¬ 
tion of quantum states of light through optical frequency 
conversion mm, or measurement-induced nonlinearities 
for quantum logic gates [4], together with linear opera¬ 
tions between optical modes to implement core process¬ 
ing functions [5]. Encoding qubits in the polarisation of 
photons has been particularly appealing for the ability 
to implement arbitrary linear operations on the two po¬ 
larisation modes using a series of wave plates [6]. For 
path encoding the same operations can be mapped to 
a sequence of beamsplitters and phase shifters. In fact, 
since any linear optical (LO) circuit is described by a 
unitary operator, and a specific array of basic two-mode 
operations is mathematically sufficient to implement any 
unitary operator on optical modes [7], it is theoretically 
possible to construct a single device with sufficient versa¬ 
tility to implement any possible LO operation up to the 
specified number of modes. 

Here we report the realisation of this longstanding 
goal with a six-mode device that is completely repro¬ 
grammable and universal for LO. We demonstrate the 
versatility of this universal LO processor (LPU) by ap¬ 
plying it to several quantum information protocols, in¬ 
cluding tasks that were previously not possible. We im- 
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plement heralded quantum logic gates at the heart of the 
circuit model of LO quantum computing [4] and new her¬ 
alded entangling gates that underpin the measurement- 
based model of LO quantum computing Hfln], both of 
which are the first of their kind in integrated photonics. 
We perform 100 different boson sampling [IIHIS] exper¬ 
iments and simultaneously realise new verification pro¬ 
tocols. Finally, we use multi-particle quantum interfer¬ 
ence to distinguish six-dimensional complex Hadamard 
operations, including newly discovered examples, where 
full classification remains an open mathematical prob¬ 
lem. The results presented required reconfiguration of 
this single device to implement ^1,000 experiments. 

Isolated quantum mechanical processes, including loss¬ 
less LO circuits acting on photons, preserve orthogonality 
between input and output states, and are therefore de¬ 
scribed by unitary operators. While the relevant math¬ 
ematics for parameterising unitary matrices has been 
known for at least a century [16], the theoretical formu¬ 
lation of a LO circuit in which an arbitrary unitary op¬ 
erator could be realised is more recent [7] . The full space 
of m-dimensional unitary matrices (or unitaries) can be 
parametrised as a product of ~ m‘^/2 two-dimensional 
unitary primitives, each with two free parameters. Any 
given unitary then corresponds to a unique set of param¬ 
eter values. An arbitrarily reconfigurable LO network 
would be realised by implementing these unitary primi¬ 
tives as two-mode Mach-Zehnder interferometers (MZIs) 
with two beamsplitters (or integrated directional cou¬ 
plers) and two phase shifters, such that any given LO 
network corresponds to a set of phase shift values. 
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FIG. 1. Universal linear optical processor (LPU). (a) Decomposition of a fully parametrised unitary for an m-mode circuit 
to realise any LO operation. Sub-unitaries Di consist of Mach-Zehnder interferometers Mij (MZIs) built from phase shifters 
(yellow) and beam splitters, to control photon amplitudes {aij) and phases (b) Multi-photon ensembles are generated 

via spontaneous parametric down-conversion (SPDC), comprising a BiBO crystal, dichroic mirrors (DM) and interference filter 
(IF); preceded by a pulsed Ti:sapphire laser and second harmonic generation from a BBO crystal. Photons are collected into 
polarisation maintaining fibers and delivered to the LPU via a packaged v-groove fibre array. The processor is constructed over 
six modes as a cascade of 15 Mach-Zehnder interferometers, controlled by 30 thermo-optic phase shifters, set by a digital-to- 
analgoue converter (DAG) and actively cooled by a Peltier cooling unit. Photons are then out-coupled into a 2nd packaged 
VGA and sent to 6 (or 12 with fibre splitters for single-mode photon number resolving capability) single photon avalanche 
diodes (SPADs) and counted using a 12-channel time-correlated single-photon counting module (TGSPG). See Appendix for 
further details. 


Such operation can be understood, with reference to 
the schematic shown in Fig. [^a), as a sequence of m 
sub-unitaries Di^ each of which is guaranteed to enable 
transformation of the state of a photon input in its top 
mode into an arbitrary superposition of all its output 
modes. The top MZI in each Di is set to retain the 
desired probability for a photon occupying the top mode, 
before imparting the desired phase on the top mode. The 
remainder of the photon undergoes a similar operation on 
the subsequent modes. Full reconfigurability is realised 
by feeding all output modes from one Di> into all but the 
uppermost input mode of the following Di>_i HZ]. 

Realising such a scheme requires sub-wavelength 
stability and high fidelity components to support both 
classical and quantum interference — possibilities 
opened up by integrated quantum photonics plSUSb] . A 
schematic of our LPU is shown in Fig. f^b). The device, 
made with planar lightwave circuit (PLC) technology 
[271 EH], comprises an array of 30 silica-on-silicon 
waveguide directional couplers with 30 electronically 


controlled thermo-optic phase shifters, to form a cas¬ 
cade of 15 MZIs across six modes (see Appendix for 
further details). Such an LPU can implement any LO 
protocol, and here we focus on several protocols at the 
forefront of quantum information science and technology. 

Quantum gates: With the addition of single pho¬ 
ton sources and measurements, and rapid feed-forward 
of classical information, both the circuit [4] and 
measurement-based isHini models of digital quantum 
computing can be efficiently implemented with LO. Ba¬ 
sic two-qubit processes are realised probabilistically with 
LO circuitry, therefore a key requirement for scalability is 
that a successful operation is heralded by the detection of 
ancillary photons, to signal that the processed photonic 
qubits are available for use in the larger architecture [29| . 

We programmed our device to implement a new com¬ 
pact four-photon scheme, suitable for measurement- 
based quantum computing, which generates the entan¬ 
gled state of two photonic qubits upon the detection of 
another two ancilla photons, schematically displayed in 
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FIG. 2. Gates for linear optical quantum computing, (a) The heralded Bell state generator receives four photons and emits two 
of them in a maximally entangled state upon detection of the remaining two. (b) Our linear optical protocol emits a Bell state 
on modes {a,b} with input and heralding modes labelled by |1). (c) Experimental data measuring photon correlations (blue), 
with error bars from Poissonian counting statistics. The overlaid ideal theoretical values show that measurements in common 
bases should be correlated while measurements in different bases should be uncorrelated (see Appendix for details), (d) The 
heralded GNOT gate is successfully implemented on the two photonic qubits, upon detection of the two ancilla photons, (e) The 
linear optical protocol realising the heralded GNOT operation on the control co,i and target to,i qubits, (f) Experimental data 
showing the computational truth table, with the ideal theoretical truth table overlaid, (g) Quantum process tomography of an 
unheralded (or post-selected) two-qubit GNOT gate can be performed with the addition of arbitrary single qubit preparation 
and measurement operations, (h) The linear optical circuit realising the post-selected GNOT gate, with MZIs (inset) allowing 
single qubit operations, (i) Experimentally determined process matrix with ideal theoretical overlaid, (j) An arbitrary single 
qubit operation can be realised with an MZI and additional phase shifters. Three consecutive MZIs allows us to perform full 
process tomography on any single qubit operation, (k) Experimental data showing the measured process matrices for the three 
Pauli operations, the Hadamard gate (H), and the tt/S phase gate (T). Experimental data is corrected for measured detector 
efficiencies. 
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Fig. ,b). For Bell states, measurements in common 
bases should be correlated while measurements in differ¬ 
ent bases should be uncorrelated. We implemented both 
types of measurement, data shown in Fig, [^c) , finding 
the mean statistical fidelity Ts = between 

all four experimental and theoretical probabil¬ 
ity distributions to be Tg = 0.966 ± 0.004 (see Methods 
for an explanation of all error analyses). We used these 
measurements to verify the entanglement of our state by 
calculating E = l/2{{ax 0 cfx) E {cFz ^ ^z))i finding a 
value of = 0.673 ± 0.031 where E > 1/2 witnesses 
entanglement m- 

Next, we reprogrammed the device to realise a quan¬ 
tum logic gate suitable for use in the circuit model of 
quantum computing [Figs. Hd, e)]. The detection of 
two ancilla photons heralds the implementation of a 
controlled-NOT (CNOT) operation [11 [311 [32] between 
two photonic qubits: the logical state of the target qubit 
is flipped (0 1) if the control qubit is in the state 1, and 

left unchanged if the control qubit is in state 0. We mea¬ 
sured the logical truth table for this operation, shown in 
Fig. if), and found its mean statistical fidelity averaged 
over all computational inputs to be = 0.927 ± 0.004 
with the ideal. 

While these are the first heralded two-qubit processes 
in integrated optics, several examples of photonic chips 
specifically fabricated to implement unheralded (and 
hence not scalable) two-qubit gates have been reported 
[2ll|33l|34]. To compare the performance our universal 
processor against devices fabricated for a specific task, we 
implemented an unheralded two-photon CNOT gate |35F 
EZ] with single qubit preparation and measurement ca¬ 
pabilities [Figs. [^g,h)] and performed full quantum pro¬ 
cess tomography[38] [Fig.j^i)]. The process fidelity was 
found to be Ep = 0.909 ±0.001, and the average gate 
fidelity Eg = 0.927 ± 0.001 (see Appendix for details), 
greater than those previously reported [2111331139]. 

Combined with two-qubit operations, a small set of 
single qubit gates [Fig. [^j)], including the Hadamard 
{H) and tt/S (T) gates are sufficient to reaiise a uni- 
versai gate set for quantum computing [40] . We impie- 
mented and performed fuii quantum process tomogra¬ 
phy for these two gates and for the three Pauii gates, as 
shown in Fig. [^k), finding an average process fideiity of 
Ep = 0.992 ± 0.008. 

Note, in muiti-photon experiments deviations from 
unit fideiity are primariiy caused by imperfections in 
the photon source such as reduced quantum interfer¬ 
ence between different pair creation events [41], and 
higher order terms in the SPDC process. To omit 
these effects and measure the performance of our LPU 
directiy we used one and two-photon ensembies to 
recover the impiemented unitary U' [42] and caicuiated 
the unitary fideiity Eu = |Tr(l/’*'.l/')/6|^ with the 

intended unitary U. For the Beii state generator, 
the heraided- and unheraided-CNOT respectiveiy we 


found Eu = 0.998±0.002,0.994±0.008 and 0.988±0.005. 

Implementing boson sampling: The reaiisation of 
a iarge-scaie quantum computer that demonstrates an 
intrinsic exponentiai advantage over ciassicai machines 
wouid be in confiict with a foundationai tenet in com¬ 
puter science: the extended Church Turing thesis (ECT). 
The ECT conjectures that aii reaiistic physicai systems 
can be efficientiy simuiated with a probabiiistic Turing 
machine, or ciassicai computer. Designed for LO, with no 
requirement for quantum iogic gates or qubit entangling 
operations, boson sampling mHis] is a quantum protocol 
that has been developed as a rapid route to challenge the 
ECT and demonstrate that quantum physics can be har¬ 
nessed to provide fundamentally new and non-classical 
computational capabilities. 

Based on the foundations of computer science, boson 
sampling is a mathematical proof (using plausible con¬ 
jectures) that a many-photon state, when acted on by a 
large LO circuit set to implement a Haar-random unitary, 
will give rise to a probability distribution that cannot be 
efficiently sampled by a classical algorithm. Quantum 
interference among the photons [43] contributes to the 
pattern of the probability distribution. The classical in¬ 
tractability arises because the probability amplitude for 
each correlated photon detection event is given by a clas¬ 
sically hard function, known as the permanent [44], of 
the submatrix that describes a particular route of pho¬ 
tons through the circuit. Experimentally, each detection 
event represents a sample drawn from that classically for¬ 
bidden probability distribution. 

Acting on three-photon ensembles we programmed our 
device to implement 100 different boson sampling rou¬ 
tines. Each circuit configuration was chosen randomly 
from the Haar measure, which we implemented via a 
direct parameterisation of our circuit parameters [45], 
the probability density functions of which are displayed 
in Eig. I^a). Eor each implementation [Eig. [^b)], we 
counted detection events for each of the 20 collision free 
ways in which three photons can exit the six output 
ports of the device, and in Eig. [^c) plot a histogram 
of fidelities with statistics based on calculations of 
matrix permanents, finding a mean statistical fidelity 
of Eg = 0.950 ± 0.020. These results show a test of the 
performance of our LPU over many circuit configurations 
randomly and unbiasedly chosen from the full space of 
all possible configurations. 

Verifying boson sampling: An open and important 
question, particularly in light of the ECT, is how to verify 
that boson sampling continues to be governed by the laws 
of quantum mechanics when experiments reach the scale 
that classical computers can no longer simulate [464149] . 
Unlike certain algorithms for digital quantum comput¬ 
ers, including Shor’s factoring algorithm [50], where the 
solution to a classically hard problem can be efficiently 
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Boson Sampling 




FIG. 3. Boson sampling and Complex Hadamard Matrices, (a) A Haar random unitary can be implemented by choosing 
beamsplitter reflectivities (or equivalently MZI phases a (inset 1)), and phase shifter values 0 (inset 2) from the probability 
density functions (pdfs) corresponding to linear optical elements in (b). (c) A histogram of experimentally determined statistical 
fidelities for 100 three-photon boson sampling experiments, with one and two photon histograms inset, (d) Experimental points 
(red) showing measured three-photon violations (ly) of the ZTL whilst scanning across the two-parameter O 2 ) matrices. 

Minima occur at the Fourier transform point {Fq) as predicted by the theoretical manifold, black lines to guide the eye. (e) After 
15 steps, dynamic updating with Bayesian model comparison determines with a confidence > 99% that six-photon detections 
(from the input state 181 , 32 )) are sampled from an indistinguishable (quantum) as opposed to a distinguishable (classical) 
distribution. (f,g) The two parameter family of 6-dimensional CHMs F^^ can be directly realised on optical modes, however 
when mapped to our general purpose circuit, changes to these phases (^ 1 ,^ 2 ) require all on-chip phase shifters to be changed. 
(h,i) Two-photon data for a single example of the newly discovered 4-parameter family G^q \ which is not a member of F^\ and 
the isolated case, with ideal theoretical data overlaid (error bars assume Poissonian counting statistics), (j) Two photon 
correlation manifolds for Fq^\ with experimentally measured probabilities (red) as colour coded in (g). 




























































































































6 


classically checked, there seems to be no analogous way 
to check that large-scale boson sampling is sampling from 
a probability distribution that arises from many-photon 
quantum interference. 

While it is likely that boson sampling is in princi¬ 
ple mathematically unverifiable, methods have been pro¬ 
posed to gather supporting or circumstantial evidence 
for the correct operation of the protocol [li[5ll|52]. The 
essence of these methods is to implement experiments 
that share basic quantum mechanical features with bo¬ 
son sampling, but where certain properties of the exper¬ 
imental output can be predicted and therefore checked. 
The zero transmission law (ZTL) [53] predicts that cor¬ 
related photon detection for most of the exponentially 
growing number of configurations is strictly suppressed 
if the circuit is set to implement the Fourier transform 
(FT) on optical modes. This is known because the struc¬ 
ture of the FT allows these matrix permanents to be ef¬ 
ficiently evaluated without explicit calculation. Because 
large scale many-photon quantum interference is at the 
core of the ZTL, it has been proposed as a certificate for 
the capability of a device to implement boson sampling 

m- 

We programmed our device to implement 16 examples 
of the Fq T^i,^ 2 ) two-parameter set of six dimensional 
matrices, including Fq (the 6 dimensional FT), which oc¬ 
curs at 6 >i, 6>2 = 7 r, 0 . Using statistics from three-photon 
ensembles, we calculated the experimental violation 
of the ZTL as z/ = Ng/N, the ratio of the number of 
predicted suppressed events Ng to the total number of 
events N; the results of which are plotted in Fig. [^d) 
alongside the theoretical manifold. The experimental 
points follow the shape of the manifold with the minimal 
violation of the ZTL Umin = 0.319 ± 0.009 occurring 
when Fq is implemented. The average ZTL violation of 
the nine points that are predicted to maximally violate 
is i^max = 0.638 ± 0.029. Crucially, this verification 
protocol is implemented in the same device and with 
the same procedure as that for the boson sampling 
experiments above. 

Six photon verification: An essential requirement of 
boson sampling is the indistinguishability among pho¬ 
tons. With our device set to implement Fe, we in¬ 
jected the six-photon state | 3 i, 32 ) and counted six- 
photon statistics with an all-fibre beam-splitter between 
each output mode and two SPADs to give probabilis¬ 
tic number-resolved photon detection over a total of 12 
SPADs. 

We used Bayesian model comparison (see Ref. 021 ) to 
update, in realtime, confidence that we are sampling from 
a (pre-calculated) quantum probability distribution (aris¬ 
ing from completely indistinguishable photons) or from 
a classical probability distribution (arising from com¬ 
pletely distinguishable photons); as shown in Fig. [^e). 
After collecting 15 six-fold coincidence events, we deter¬ 
mine a confidence of p = 0.998 that we are sampling from 
a quantum (not classical) distribution. 


( 2 ) 

Complex Hadamard operations: The FT and Fq ^ 

are examples in the more general class of complex 
Hadamard matrices (CHMs), which are related to mutu¬ 
ally unbiased bases [54| and are of fundamental interest 
in quantum information theory m- CHMs are defined 
as X unitary matrices with entries of squared abso¬ 
lute value equal to 1/N. While this definition is straight¬ 
forward, classification of these matrices is far from trivial 
and is concerned with identifying CHMs that are inequiv¬ 
alent up to pre and post multiplication with permutation 
matrices and diagonal unitaries [56]. In the N = {2, 3, 5} 
case, all CHMs are equivalent to the respective FT ma¬ 
trix, while for N = 4 there exists a one parameter equiv¬ 
alence class. Whilst a full classification of N = 6 CHMs 
is unknown, it is currently conjectured that the set con¬ 
sists of an isolated matrix which does not belong to 
any family 03, and a newly discovered four-parameter 
generic family [58] . 

In LO experimental implementations, discrimination 
among CHMs can be accomplished via the observation of 
characteristic patterns of photonic quantum interference 
[59f[62] . Up until now, these observations have been too 
experimentally challenging for the six-dimensional case. 
Here, we reconstruct correlation manifolds of two-photon 
detection events by scanning over the Fq ^ matrices, dis¬ 
playing four (out of the 15 sets collected) in Fig.j^j). We 
find a mean statistical fidelity of Fg = 0.979 ± 0.007. 

We then implemented an instance of G^q^ (that is 
not contained in Fq^^) and Sq^\ and observed their 
predicted characteristic two-photon quantum inter¬ 
ference patterns displayed in Figs. ih, i) respectively, 
finding statistical fidelities to be Fg = 0.986 ±0.001 and 
Fg = 0.998 ±0.001. The intractability of calculating the 
permanents of certain CHMs is an interesting research 
line, as is the possibility of searching for new CHMs 
using photonic statistics. 

Discussion 

Photonic approaches to quantum information science 
and technology promise new scientific discoveries and 
new applications. Linear optical circuits he at the heart 
of all of these protocols, and a single LPU device with the 
ability to arbitrarily ‘dial-up’ such operations promises 
to replace a multitude of existing and future prototype 
systems. In the near term, combining LPUs with exist¬ 
ing higher efficiency sources and detectors will expand 
their capabilities, while in the mid-term, integration of 
these components [63] [64] with larger low-loss circuits 
[65] will open up new avenues of research and application. 

Methods: Circuit implementation analysis 

LO protocols for quantum information processing are 
typically conceived and explained as a network of 
optical modes connected with beamsplitters and phase 
shifters. Provided the total number of modes in such a 
protocol is not greater than that of an m-mode LPU, 
the network can be mapped to the device via its unitary 


7 


matrix description, or (due to symmetry and structure) 
a direct implementation can be recognised. In the 
following sections we explain the circuits implemented, 
and note that both compiling methods are used. While 
the generality of this LPU makes it well suited to 
implementing random unitaries, it is not optimised for 
the highly structured unitaries that typically implement 
qubit operations. Nevertheless, this system demon¬ 
strates fidelities greater than those reported for optical 
chips that have been custom made for LO quantum 
processing tasks. Throughout this manuscript two error 
analysis methods are used: those for individual fidelities 
T calculated by propagating Poissonian count rate 
errors, and those for mean fidelities T calculated as lajr. 
The bar symbol will used to denote which method is used. 

Heralding entangled states: Our compact Bell state 
generator (BSG) scheme shown in Fig. nib), accepts 
four photons into modes {1,3,4, 5}. Detection in ancilla 
modes {3,4} occurs with probability p = 2/27 and her- 
aids the state |i>+) = (|li, 02 , l 5 , 06 ) + |0i, 12 , 05 , le))/V2 
(where |nj) represents n photons in the i**’ mode) over 
the remaining modes with unit probability in the ideal 
case. Operating the BSG with the non-ideal state from 
our SPDG source requires the Bell state to be detected 
to remove the unwanted photon terms (see Appendix for 
details). 

After implementing the BSG, measurements on both 
qubits in the computational basis G are straightfor¬ 
ward and arbitrary single qubit measurements on the 6o,i 
modes (Jz^cFn are possible with an unused MZI. We were 
also able to measure both qubits in the conjugate basis 
G (Jx by adding a fibre-looped Sagnac interferometer 
on the ao,i modes to send this qubit back through the 
chip to meet an MZI at the input facet (see Appendix 
for details). 

We found the statistical fidelity for the common basis 
measurements (Tz^cFz and (Tx^cFx to be Tg = 0.891 ±0.015 
and 0.979 ± 0.003 respectively, and the uncommon bases 
(Jz'^cFx and (Jz'^cFy to be 0.999±0.002 and 1.000±0.002. 

Heralded quantum logic gate: The heralded CNOT 
gate shown in Fig.j^e), is a four-photon scheme [32] that 
has been adapted for implementation with the six-mode 
LPU, plus a seventh auxiliary coupling mode (mode 0) 
also contained on our chip. The function of an 8th lossy 
balancing mode is implemented in post-processing as an 
equivalent reduction in detection efficiency on the lower 
ancilla mode (see Appendix for details). 

Control and target qubits co,i and to,i are imple¬ 
mented via the modes {0,2,3,4} respectively, and 


heralding modes via {1,5}. In the ideal case, the 
detection of the two ancilla photons, which occurs with 
probability p = 1/16, guarantees the CNOT logic has 
been performed on the two photonic qubits. Once again, 
the non-ideal input state requires a computational 
detection to remove unwanted terms. 

Full quantum process tomography for 1 and 2 
qubits: The unheralded (or post-selected) CNOT gate 
shown in Fig. [^h) is a two-photon scheme that requires 
the detection of the two photonic qubits in a computa¬ 
tional state to signal the success of the gate with proba¬ 
bility p = 1/9 [351137] . As the scheme requires only two 
photons (which our source produces with high fidelity) it 
allows for a less obfuscated test of our device for logical 
operations. 

Computational states co,i and to,i are implemented 
via modes {2,3,4,5}, and MZIs (shown in orange in 
Fig. I^h)) allow arbitrary product state preparation and 
measurement of both qubits, so truth tables may be taken 
in any bases. We perform full quantum process tomogra¬ 
phy using an over complete set of 324 truth tables from 
all possible combinations of Pauli eigenvector input states 
and Pauli basis measurements (see Appendix for details), 
finding a mean statistical fidelity between all experimen¬ 
tally reconstructed truth tables and ideal values to be 
fs = 0.986 ±0.001. 

We performed single qubit process tomography on 
the Hadamard (i7), n/S (T) and (Jx^cFy^cTz gates, 
finding process fidelities of 0.995 ± 0.001, 0.978 ± 0.001, 
0.997 ± 0.001, 0.993 ± 0.001, 0.994 ± 0.001 respectively, 
demonstrating very high fidelities for operations requir¬ 
ing only one-photon interference. 

Haar random unitaries: Along with three-photon 
ensembles injected into modes {1, 2, 3}, as a benchmark¬ 
ing exercise we injected the one and two photon states 
|li) and 11 1 , 12 ) yielding fidelities of 0.999 ± 0.001 and 
0.990 ± 0007 respectively. For three of our three-photon 
boson sampling unitaries (those that gave the highest, 
lowest, and mode fidelity) we recalculated the statistical 
fidelities based on the recovered unitary U' . This 
increased the mean statistical fidelity and reduced range 
jpmax _ jpmin q ^^2 to 0.023; implying dialling 

errors are non-uniform across the space of unitaries. 

Compiling : Figure ^g) shows a circuit customised 
to implement the Fq ^ family of CHMs, with the defining 
parameters ^ 1,^2 explicitly realised. However due to the 
particular decomposition of our device, compiling F^^ 
onto our LPU requires unique phase shift values across 
the entire device for each implementation ^ 1 , ^2 ^ [0, 27r). 
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Appendix A: UNIVERSAL LINEAR OPTICAL 
PROCESSOR 

1. Photon Source 

As shown in Fig. [^b), SOSnm laser light from a 140fs 
pulsed Titanium: Sapphire laser was attenuated via a 
half wave plate (HWP) and polarising beam splitter 
(PBS), before being frequency doubled via a type-I Bar¬ 
ium borate (BBO) nonlinear crystal. The subsequent 
404nm light was reflected from dichroic mirrors (DM), 
filtering the 808nm light into a beam dump (BD), and 
focused onto a type-I Bismuth Borate (BiBO) nonlinear 
crystal to generate multiple pairs of photons through 
spontaneous parametric down conversion (SPDC). After 
passing through an interference filter (IF) photons are 
collected into four polarisation maintaining fibres (PMF) 
and delivered to the LPU. 


2. Fabrication 

The silica-based planar waveguide circuit is fabricated 
on a silicon wafer using flame hydrolysis deposition, pho¬ 
tolithography and reactive ion etching. The germanium 
(Ge 02 ) doped silica core has a cross sectional dimen¬ 
sion of 3.5/im X 3.5/im, and is surrounded by a silica 
cladding yielding a core-cladding refractive index differ¬ 
ence of 0.45%. The input and output waveguides are 
connected to polarization-maintaining fibre arrays by us¬ 
ing ultraviolet curing epoxy adhesive with in- and out- 
coupling losses of approximately 0.4dB/facet (9%). 

The circuit consists of 30 tunable phase shifters, each 
phase of which is controlled by thin-film thermo-optic 
heaters fabricated on top of the circuit. The electric 
power required for a 27r phase shift is ^ 0.81F. Direc¬ 
tional couplers are designed to have a length of 500/im 
and a waveguide gap of 2 /im, with each directional cou¬ 
plers loss estimated to be less than 0.02dB (0.4%). The 
mean insertion loss (fibre to fibre) when averaged over all 
modes is 2.4dB (42%). 

3. Electrical Control 

We developed a custom electronic circuit to deliver 
power to each of the on-chip heaters independently. This 
circuit consisted of four sets of the following: one micro¬ 
controller for communication, arithmetic, and measure¬ 
ment; eight 12 -bit digital-to-analogue converters; eight 
high-power non-inverting drive amplifiers; eight current- 
measurement amplifiers; and eight 12 -bit analogue-to- 
digital converters for reading out the current sourced to 
each heater. Each of the 32 drive ports delivered up to 
20 V, with 4.9 mV resolution, and could source up to 
100 mA of current, with a measurement resolution of 
24 /iA. 


4. Circuit Characterisation 

Each thermo-optic phase shifter is described by a non¬ 
linear phase-voltage relationship 

$(y) = a ++ 7y^ (Al) 

where 4> is the resulting phase shift due to an applied 
voltage V and a, (3 and 7 are real numbers, which in gen¬ 
eral vary for each heater due to imperfections in the fab¬ 
rication process. To determine 4>(V) we inject heralded 
single photons into the fundamental transverse-magnetic 
(TM) mode of an isolated MZI (<ai,j) or phase shifter 
( 0 ij) then sweep the voltage 1.8 ^ 10 V producing an 
interference fringe given by 

C = A-5.cos^(V), (A 2 ) 

where C is the count rate in the cross mode and A and 
5, along with a, (3 and 7 are free parameters to be de¬ 
termined via a least squares fitting procedure. 

The triangular arrangement of our device allows each 
phase shifter to be isolated and characterised indepen¬ 
dently of any other. We first characterise the diagonal 
0^17 (defined in Eigurej^) by: 

1 . injecting single photons into mode 1 

2 . measure interference fringe on 0^17 

3. set Qfiq = 27r 

4. repeat for ai^j for j = [2, 5] 

completing the characterisation of ai^j. We repeat this 
procedure for aij^ i = [2,5] by injecting photons into 
mode i and setting all = tt, thus completing the 

characterisation of all aij. We then characterise (j)i^^ by: 

1. injecting single photons into mode 5 

2. set 0^75 = TT for j = [3, 5] 

3. set 0^74 = TT for j = [1,4] 

4. set 0 ^ 2 , 5 , 0 ^ 1,5 = 7r/2 

5. measure interference fringe on ^ 1^5 

6 . repeat for ^ 75 , with ai+ 1,5 = 7r/2 and i = [2,4] 

completing the characterisation of ^ 75 . We then repeat 
this procedure for j = [ 2 ,4] by injecting photons into 
mode j and setting = tt, thus completing 

the characterisation of all and therefore all phase 
shifters. 
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5. Cross Talk 

Electrical cross talk was observed due to certain 
heaters sharing a common ground. This is corrected for 
in pre-processing by fully characterising all shared resis¬ 
tances and calculating the set voltages {V-} which yields 
the actual voltages {Vi} required for a desired phase shift. 
Once corrected for no evidence of cross talk, thermal or 
otherwise, is observed. 

A final correction is applied due to an observed reduc¬ 
tion in quantum interference at certain points within the 
chip caused by polarisation rotation. We correct for this 
by inserting polarising beam splitters at the output of 
the chip filtering out undesired polarisation states. 


where K is the number of outcomes in the measurement 
of pij and /3 = 0.1 is the hedging parameter. The min¬ 
imisation is performed by phrasing it as a Semidefinite 
Program (SDP) and solved using the CSDP solver j67] . 
The average gate fidelity, J'g, defined as 

was calculated from the 
process fidelity via the relation Tg = {dTp + l)/(d + 1) 
All errors were calculated by a bootstrapping method, 
where the raw experimental counts were resampled 100 
times with Poissonian noise and the statistics of the re¬ 
sulting distributions analysed. 


Appendix C: QUANTUM LOGIC GATES 


^ 1. Unheralded CNOT 

6. Phase accuracy 


To determine the accuracy in setting a phase we run 
a randomized benchmarking type experiment by inject¬ 
ing heralded single photons into mode 1, then setting 
100 different vectors dij = 0 ^ 1 , 2 , <^ 1 , 3 , <^ 1 , 4 , <^ 1 , 5)5 

chosen from the Haar measure, to give a mean fidelity 
yrexp between our experimental and theoretical distribu¬ 
tions. We then run a Montecarlo simulation of the ex¬ 
periment by applying Gaussian noise (5^ to each phase 
shifter and fitting this fidelity as a function of this 
noise, giving We then solve 

finding (5^ = 0.035 rad. This value lumps together error 
caused by the circuit, the control electronics and the cal¬ 
ibration procedure, so is a good metric for our effective 
phase accuracy. 


Appendix B: QUANTUM PROCESS 
TOMOGRAPHY 


All quantum process tomography was performed as fol¬ 
lows. Measurement effects were of the form \{)n) {{>n \ ^ 
\{’m) {{’m\ and preparations were of the form \{>n) ^ \{’m) 
where \^l;) G {|0) , |1) , |+) , |-) , |-hi) , |-^)}. An estimate 
of the Choi-Jamiolkowski state corresponding to the 
process was found using maximum-likelihood estimation 
(MLE) by minimisation of an approximation to the neg¬ 
ative loglikelihood function: 

est • / ~ dTl [psEij])‘^ I 

—1 

subject to: pg > 0, Tr^ (pg) = 1/d 


where pij = riij/Nij are the experimental frequencies 
and Eij = {piY i labelling the preparation and 

j the measurement effect. To avoid problems associated 
with zero probabilities in MLE we use hedged MLE [66] , 
where the experimental probabilities are adjusted to 


Pij = 


riij -h P 

Nij + KP 


The two photon, unheralded CNOT gate is imple¬ 
mented using modes {2,3} and {4,5} as dual-rail qubits, 
inputting photons in modes {2,4} and post-selecting on 
coincidence events in modes {{2,4},{2,5},{3,4},{3,5}}. 
To perform the gate we use the following chip phases: 


/(O.OOO, 1.231) 


R = 


(012, Q;i2) 

(0.000, q;22) 


V 


(0.000,3.142) (014, au) 

(023,1.231) (0.000,1.571) 

(0.000,3.142) (3.142,1.571) 

(0.000, q;44) 


(0.000,3.142)\ 
(3.142,3.142) 
(0.000,1.231) 
(045,3.142) 
(0.000,3.142)/ 


Where Rij = {(j)ij,aij). { 0 ^ 22 , ^ 23 } and { 0 ^ 44 , ^ 45 } can 
be set to prepare arbitrary pure states of qubits 1 and 
2 respectively, and { 0 i 2 ,<ai 2 } and { 0 i 4 ,(ai 4 } can be set 
to perform arbitrary projective measurements on qubits 
1 and 2 respectively. To prepare Pauli eigenstates and 
measure Pauli observables, the following phases are set: 


State 

a 

0 

Observable (j) 

a 

|0) 

TT 

0 

CTz 

0 

TT 

ll) 

0 

0 


0 

TT 

2 

1+) 

TT 

2 

0 

(Jy 

TT 

2 

TT 

2 

I-) 

TT 

2 

TT 




1+0 

TT 

2 

TT 

2 




1-0 

TT 

37r 




2 

2 




2. 

Heralded CNOT 



The unitary operation (up to input and output phases) 
required on our six mode circuit to perform the heralded 
CNOT gate is given by 


/ 0.476 

-0.622 

-0.440 

0.440 

0 

0 \ 

-0.622 

-0.476 

0 

0 

0.622 

0 

-0.383 

0 

0.293 

0.707 

-0.383 

-0.348 

0.383 

0 

0.707 

0.293 

0.383 

0.348 

0 

0.622 

-0.440 

0.440 

0.476 

0 

\ 0.306 

0 

0.166 

-0.166 

0.306 

-0.870/ 
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Corresponding to a set of phases 


/(O.OOO, 0.992) 


R = 


(4.712,1.571) 
(0.000,1.571) 


(4.544,4.957) 

(1.571,0.992) 

(0.000,1.571) 


V 


(5.375,1.792) (3.816,0.000)\ 

(2.188,0.000) (4.712,0.000) 

(5.498,0.000) (1.571,2.226) 

(0.000,3.142) (3.142,0.000) 

( 0 . 000 , 0 . 000 )/ 


The truth table of the gate was experimentally tested by 
inputting four photons, two in the ancilla modes and two 
in the computational modes. Our photon source is based 
upon two independent SPDC events which produce the 
state ^^ch that in the four photon subspace 

we have |'04) ^ + |2200) + |0022). Although we do 

not herald the ideal input the other terms should 

not produce output states in the computational subspace 
and so could be removed in measurement for this proof- 
of-principle demonstration. 


3. Bell state generator 


The circuit used for Bell state generation was given by 
the unitary matrix: 


U = 


/0.707 0.707 0 

0.408 -0.408 -0.577 

0.408 -0.408 0.289+ 0.289i 
0.408 -0.408 0.289-0.289i 
0 0 0.333 - 0.471i 

V 0 0 0 


0 

0.577 

-0.289 + 0.289i 
-0.289 - 0.289i 
0.333 - 0.471i 
0 


0 

0 

-0.408i 

0.408i 

0.236 - 0.333i 
0.707 


0 \ 

0 

-0.408i 

0.408i 

0.236 - 0.333i 
0.707 / 


Corresponding to a set of phases 


/(O.OOO, 1.571) 


R = 


(0.000,1.231) 

(0.000,3.142) 


(0.000,1.571) 

(0.000,1.571) 

(0.000,3.142) 


V 


(0.000,3.141) 

(1.571,1.231) 

(0.000,3.142) 

(0.000,3.142) 


(0.000,3.142)\ 

(0.000,3.142) 

(0.000,1.571) 

(0.000,3.142) 

(0.000,3.142)/ 


To obtain statistics for the ax measurement on qubit 1, 
we create a fibre looped Sagnac interferometer between 
the output of modes uq and ai (see Fig. 2b) injecting 
the photon back into the chip with a stable phase. Due 
to the configuration of the circuit, the photon sees a 1/3 
reflectivity MZI between modes 2 and 3, before a 1/2 
reflectivity MZI between modes 1 and 2. This circuit 
implements a single qubit POVM with three elements, 
, E _, El} where 

E+ = i|l)(l| + i|0)(0| + ^(|l)(0| + |0Xl|) 
i;_ = i|i)(i| + i|o)(o|-^(|i)(o| + |o)(i|) 

Ei = l | 1 )( 1 | 

corresponding to detection in modes {1,2,3} respec¬ 
tively. Since we cannot detect in mode 3 (this is one 
of the inputs), we postselect on events where the photon 
is detected in modes 1 and 2, which should occur with 
probabilities given by 


P- 


T^{E+p) 

Tv{E+p) + Tr{E_p) 
Tr {E-p) 

Tv{E+p) + Tr{E_p) 


which leads to an expectation value 

p+-p- = (a,) = iTr(p|i)(i|)+Tr(p|0)(0|) 

such that, with knowledge of the data from the mea¬ 
surement, we are able to calculate the expectation of 

Since the experiment would require both input and 
detection in mode 1, we implement it in two stages, first 
inputting in mode 1 and detecting four-fold coincidences 
in modes {2,3,4,5} and {2,3,4,6}, then inputting in mode 
2 and detecting four-fold coincidences in modes {1,3,4,5} 
and {1,3,4,6} (note that by inputting in mode 2, the state 
would pick up a relative tt phase shift so this is offset by 
setting 035 = ^). The two data sets are then normalised 
with respect to each other by the ratio of total counts in 
each experiment. 

As in the case of the heralded CNOT gate, it was pos¬ 
sible to remove events caused by the non-ideal input from 
the measurement statistics. 


Appendix D: COMPLEX HADAMARD 
MATRICES 


1. Zero-Transmission Law 

The most general form of the ZTL applies to any peri¬ 
odic input state to the FT matrix [68]. For an n photon, 
m-periodic initial state r input into the FT, final states 
s are suppressed when the sum of the mode assignment 
list d{^ (i.e. the list of positions of the photons in the 
output state) multiplied by m is not divisible by n 

moderndi(s),7^0 Prob(s) =0 

i 

For our implementation, we input the periodic state 
r = {1, 0,1, 0,1, 0} into the FT, and measure three-fold 
coincidences in the collision-free subspace. In this sub¬ 
space, 12 out of the total of 20 outputs are suppressed 
and u is calculated as the ratio of events in these 12 out¬ 
puts to the total events recorded. 

2 . Si°'’ &L 

The isolated CHM whose output distribution is la¬ 
belled in Fig. 3(i) is given by 
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1 
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1 
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w 

to 
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1 

1 

w 

1 

w‘^ 

w‘^ 

w 

^/6 

1 

w 

w‘^ 

1 

w 

2 


1 

to 

w‘^ 

w 

1 

w 



w‘^ 

w 

w‘^ 

w 

1 


Where w = exp(27ri/3). 


The CHM from the generic four-parameter whose dis¬ 
tribution is labelled in Fig. 3(h), is given by 

0.408 0.408 0.408 0.408 0.408 

0.085+ 0.399i -0.373 - 0.165i -0.389 - 0.124i 0.137 - 0.385i 0.303 + 0.274i 

-0.121-0.39i 0.407 +0.028i -0.196 + 0.358i -0.252 - 0.32H -0.247 + 0.325i 

0.333 - 0.236i -0.408+ 0.021i 0.076 + 0.401i 0.349 +0.212i -0.093 - 0.398i 

0.269+ 0.308i 0.37 +0.173i -0.198 - 0.357i -0.293 + 0.285i -0.019 - 0.408i 

0.4-0.081i -0.404 - 0.057Z 0.298 - 0.279i -0.35 + 0.21z -0.352 + 0.206i 


/ 0.408 
0.408 -' 
qG) _ 0.408 - 

6 0.408 -' 

0.408 -' 
V 0.408 



